arXiv:hep-ph/9706421vl 19 Jun 1997 


W UE/ITP-97-016 

|hep-ph/9706421 


Couplings of heavy hadrons with soft pions 
from QCD sum rules 


A. G. Grozin* 

Institut fur Physik, Johannes-Gutenberg-Universitat, 
Staudinger Weg 7, D-55099 Mainz, Germany 
O. I. Yakovlev* 

Institut fiir Theoretische Physik II, Bayrische-Maximilians-Universitat, 
Am Hubland, D-97074 Wurzburg, Germany 

Abstract 

We estimate the couplings in the Heavy Hadron Chiral Theory 
(HHCT) lagrangian from the QCD sum rules in an external axial field. 

We take into account the perturbative correction to the meson cor¬ 
relator in the infinite mass limit. With the perturbative correction 
and three successive power corrections, the meson correlator in an ax¬ 
ial field becomes one of the best known correlators. In spite of this, 
the corresponding sum rule is not very stable. It yields the result 
giF^/(380MeV)^ = 0.1-h 0.2, where = /^m/4 = (380MeV)^ is the 
central value of the heavy meson decay constant with the perturbative 
correction | [l^ . This result is surprisingly low as compared with the 
constituent quark model estimate gi = 0.75. The sum rules for ( 72,3 fol¬ 
lowing from nondiagonal S — S and diagonal A — S baryon correlators 
in an external axial field suggest 52,3 = 0.4-h 0.7, while diagonal S — S 
and nondiagonal A — S baryon sum rules have too large uncertainties. 


* On leave from Budker Institute of Nuclear Physics, Novosibirsk 630090, 
Russia. 





1 Introduction 


It is well known that the QCD lagrangian with ni massless flavours has the 
SU{ni)L X SU{ni)j{ symmetry spontaneously broken to SU{ni)y giving the 
(n^^ —l)-plet of pseudoscalar massless Goldstone mesons (pions) vrj (vrj = 0). 
Their interactions at low momenta are described by the chiral lagrangian 
(see e. g. Q) 

2?7r 

= + S = exp—, (1) 

^ JTV 

where the pion constant /^r ~ 132MeV is defined by 

<0\jj^,W> = 'i-UejPfj., j]^ = 9 ^ 7 ^ 75 ^* ^ 

(e*- is the pion flavour wave function), and dots mean terms with more 
derivatives. Light quark masses can be included perturbatively, and lead 
to extra terms in SU{ni)L X SU{ni)ii transformations act as S —> 

LTiR^. Let’s define ^ = expz7r//,r) 5] = it transforms as ^ = 

where C/ is a SU{ni) matrix depending on 7r(x). The vector = 
+ and the axial vector —transform as 

—> U{Vfj, + dfj,)U~^, Ua^U~^. There is a freedom in transformation 

laws of matter fields such as because it is always possible to multiply 
them by a matrix depending on tt. The only requirement is the correct 
transformation with respect to SU{ni)v {L = R). It is convenient to choose 
V' —> Then the covariant derivative = 5^ + u^ transforms as —> 

U Covariant derivatives of tensors with more flavour indices are defined 
similarly. 

Hadrons with a heavy quark are now successively investigated in the 
framework of the Heavy Quark Effective Theory (HQET) (see |p for review 
and references). To the leading order in 1/m, the heavy quark spin does 
not interact and can be rotated or switched off at all (spin-flavour and 
superflavour symmetry). The Qq mesons with a spinless heavy quark form 
the ^ n^-plet The Qqq baryons can have = 0"*“ or I"*" giving the 
scalar flavour-antisymmetric npni — l)/2-plet and the vector flavour- 
symmetric ni{ni + l)/2“plet Switching the heavy quark spin on gives 
degenerate 0~ and 1“ B and B* mesons, | A baryons, and degenerate | 
and S and S* baryons. 

Interaction of these ground-state heavy hadrons with soft pions is de¬ 
scribed by the Heavy Hadron Chiral Theory Q (see I, I for review and 
references). We start from the HHCT lagrangian with the heavy quark spin 
switched off |^: 

L = L^ + PiiDoP^ + A*jiDoA^^ + S 7 • {Wq - A)S (2) 

• a) X + 2^3 (A^.a® • + E*, • a®A^’") , 
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where A is the S-A mass difference. The possibility of consideration of 
the SAtt interaction in HHCT relies on the fact that this difference is small 
compared to the chiral symmetry breaking scale though formally both of 
them are of the order of the characteristic hadron mass scale. The matrix 
elements of the axial current between heavy hadrons are easily obtained 
using PCAC; 


<M'\fj\M> = giujjj5u\ 

<S'|J!.|S> = 2ig2e% x e*^ <A|ji.|S> = 2536*^6 


( 3 ) 


where e*-^, are the M, A, S wave functions, and the nonrelativistic 
normalization of the states and wave functions is assumed. If we switch the 
heavy quark spin on, we obtain the usual HHCT lagrangian |^. 

The HHCT couplings gi should be in principle calculable in the under¬ 
lying theory—HQET, but this is a difficult nonperturbative problem. Some 
experimental information is available on gi and g^ ■ If we neglect I /rUc correc¬ 
tions, then r(T)*''' —> = gfp^/(G tt f^), and similarly for (with 

the extra 1/2). The experimental upper limit < O.lSlMeV com¬ 
bined with the branching ratio —> D^7r~^) = 68% |^] gives gi < 0.68. 

A combined analysis of D* pionic and radiative decays |Q] gives gi ~ 0.4- 
0.7. The recent CLEO measurement Q of the S* —> Acvr decays gives jK] 
g '3 = 0.7 ± 0.1. The decay S* —> Scvr is kinematically forbidden; therefore, 
the direct measurement of g 2 is not possible. 

In the constituent quark model, gi is the axial charge g of the constituent 
light quark in the heavy meson. Moreover, following the folklore definition 
“constituent quark is B meson minus b quark”, this is the most clear way 
to define g of the constituent quark. The baryonic couplings 52,3 are also 
equal to g in this model. The most naive estimate is g ~ 1; the nucleon 
axial charge is gA = ^g in the constituent model, and in order to obtain 
gA = 1-25 we should assume g = 0.75. 

Sum rules |11| were successfully used to solve many nonperturbative 
problems in QCD and HQET. The currents with the quantum numbers of 
the ground state mesons and baryons with the heavy quark spin switched 
off are 

Jm = fh = ( 4 ) 

where Q is the spinless static quark field, C is the charge conjugation ma¬ 
trix, means q transposed, (ij) and [ij] mean symmetrization and anti- 
symmetrization. There are also currents ^rid ^*^2 with the additional 
7 o. Correlator of the mesonic currents was investigated in [^, |T^, Q, and 
of the baryonic ones—in |15, 16|. The sum rules results are in a qualitative 
agreement with the constituent quark model: the massless quark propaga¬ 
tor plus the quark condensate contribution simulate the constituent quark 
propagation well enough. 
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The sum rules method was generalized to the case of a constant external 
field for calculation of such static characteristics of hadrons as the magnetic 
moments |17]. Sum rules in an external axial field were used |18| for calcu¬ 
lation of QA of light baryons. In the present work we use HQET sum rules 
in an axial field to estimate 

We introduce the external axial field = 0) by adding the term 




( 5 ) 


to the lagrangian. We are going to calculate correlators of the currents (^) 
up to the terms linear in A (these terms are denoted by the subscript ^4). 
The light quark propagator in the gauge Xfj_Af^{x) = 0 gets the contribution 
(whose first term was found in ||l^) 


iA] ■ X 

0) = 27r2 


^75 


4x2 


+ • 


( 6 ) 


The term in Sa vanishes after the vacuum averaging; we are not going 
to calculate gluonic contributions beyond G^ and hence may omit this term. 
The axial field induces the quark condensates 




m\f; 

36 


+ 


^5x^y4 


]-2A) 




( 7 ) 


n 


<Q"^°9G^^qjbi3>A = (^j>777 - 


where mf ss 0.2GeV^ |j^, ^ is defined by <0|d(7G^i/7^u|7r’''> = 


= 3 is the number of colours, Gf = (A^^ — 1)/(2AI) (the sum rule con¬ 
sidered in 120] yields the relation ml ~ ru-Q/d). This formula has been first 
obtained in j^j, with the wrong sign of mf terms (this sign was changed 
in |21|, where a flavour-singlet external field was considered). We assume 
p ■ A = 0 where p —> 0 is the momentum of the field A; in general these 
formulae should contain Aj_ = A — [A - p/p^)p. 

The D*Df coupling was estimated ||^ from the sum rules in an external 
axial field at a finite rUc with the accuracy up to m\ term. This last term 
is incorrect; however, it is rather small. The large rUc limit was analyzed 
in @.01. We think that the treatment of the phenomenological side of the 
sum rule in is not quite convincing (see Sec. ^). In our previous work j^j, 
we used an incorrect formula for the induced condensates, which leads to 
the wrong sign of the term. After correcting this, the result agrees 


with [^, ^j. In [^1, the two-current product sandwiched between the 

vacuum and a soft pion was considered; this is equivalent to the two-point 
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sum rules in an external axial field. Both finite and infinite mass case are 
discussed in this work, with the extra power correction mQ<qq> included 
(it is rather small). Light-cone sum rules for both finite and infinite quark 
mass were considered in [^] (see also |^ 6 |). In this approach, the leading 
term of the correlator and power corrections are determined by the leading- 
and subleading-twist pion wave functions at x = 1/2. These values were 
determined from QCD sum rules |2£], and are not so accuratly known as 
the parameters of the external-field sum rules, /^r and <qq>. In the soft 
pion limit, the formulae of are reproduced (without the highest power 
correction). The 0*071 coupling (at both finite and infinite quark mass) 
was also considered in the framework of 3-point QCD sum rules, double 
Borel and double moments |28|. It is difficult to compare this approach 
with the one considered here, because there correlators are determined by a 
completely different set of contributions (with even dimensions). 

Perturbative corrections to sum rules were not discussed in the above 
works. Here we present the next-to-leading perturbative term in the infinite 
mass limit. With this term taken into account, the correlator used in the 
gi sum rule is one of the best known correlators: in addition to the leading 
term, one perturbative correction and three successive power corrections 
are known. In spite of this, the sum rule is not very stable (Sec. We 
also discuss sum rules for the baryonic couplings < 72,3 (Sec. p, following ||^ 
with some corrections. The results for diagonal A-S and nondiagonal S-E 
correlators suggest < 72,3 = 0.4 0.7. The sum rules based on nondiagonal 

A-S and diagonal S-S correlators have strong Borel parameter dependence 
and therefore yield no definite predictions. 


2 Mesons 

The correlator of the meson currents (^) has the A-term 

i<Tf{x)jjiO)>A = (8) 

that depends only on A* and not A*q. The correlator possesses the usual 
dispersion representation at any A®-^. The meson contribution at A®-^ = 0 
is p{u!) = F‘^6{io — e) where e is the meson energy and < 0 |j®|M> = iFu'' 
(the usual meson constant is /m = ‘i^Fj^fm). Switching on the external 
field produces the energy shift e ^ e — <71757 • A (s„, = 5757 • n = is 
the meson spin projection). Therefore the meson contribution is PAi^) = 
giF'^6'{u!—e)+c6{LO — s) where the second term originates from the change of 
F‘^. Besides that there is a smooth continuum contribution Thus 

we obtain 

n4H = /iiiL + _^ + nr‘H- (9) 

[(jj — ey Lo — e 
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In other words, the lowest meson’s contribution in both channels (Fig. |^) 
gives the double pole at a; = e; mixed lowest-higher state contributions 
(Fig. give a single pole at a; = e plus a term with a spectral density in 
the continuum region after the partial fraction decomposition; higher states’ 
contributions (Fig. |I|c) have a spectral density in the continuum region only. 





Figure 1: Physical states’ contributions to the correlator. 


We can also calculate the correlator using OPE. Gluons don’t interact 
with the heavy quark in the fixed point gauge x^A^{x) = 0. The leading 
perturbative diagram (Fig. ^a) vanishes. In fact, all diagrams with non¬ 
cut light quark line, including gluon condensate contributions (Fig. and 
higher loop corrections (Fig. |c) , vanish. The reason is that we can an- 
ticommute 75 out, and obtain a correlator of two currents (one of which 
contains 75 ) in a constant spatial vector field. By virtue of the Ward iden¬ 
tity, it is equal to the derivative of the correlator without external field with 
respect to the spatial momentum flowing along the light quark line. This 
is minus derivative with respect to the spatial momentum flowing along the 
heavy quark line. But the HQET propagator does not depend on the spa¬ 
tial momentum, and hence the derivative vanishes. We have checked by a 


direct calculation, using the technique of |^, that the two-loop correction 
(Fig. 1^) vanishes in the case of naive anticommuting 75 (which is the correct 
one to use in the weak current). 

The diagram with cut line (Fig. |^) gives 


nA(t) = 


f 

J 7T 


' i<qq>t 5 2,2 iml<qq>t^' 
3/2 +36"^' 48/2 


( 10 ) 


Thus the appearance of qa of the constituent quark is entirely caused by 
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Figure 2; Correlator of meson currents. 
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interaction with the quark condensate. The first three 
with the fourth one is taken from p^ ]. 

corresponds to the spectral density 


terms in ( 0 ) agree 
The correlator (K) 


PA{i^) 



<qq> 


5'{u:) 




mQ<qq> 

48/2 



( 11 ) 


We have calculated the one-loop correction to the leading dimension 
contribution. Only the one-particle-irreducible diagram (Fig. ^e) has to be 
calculated, all the other contributions can be taken from The correlator 
of bare currents in d = 4 — 2e dimensions is 


n" 


w = 




1 + 


Cpgi{-2u;)-^^ {d - 2){d - 7) 


(47r)'^/2 


2(d-3) 


r(-e)r(l + 2e) 


( 12 ) 


After renormalization of the heavy-light HQET currents, it gives the finite 
result 


11 ^( 0 ;) 


Pa{uj) 


4ll! 
f2 r 

J 71 

4 


Qlc 


l + <^i-^(- 61 og 


—2lo 


P 


+ 5 


6{u}) — 6 Cf 


^ r 

47r uj 


(13) 


In order to obtain the coordinate-space result at t = —ir, one has to cal¬ 
culate the exact d-dimensional spectral density, take its Laplace transform, 
and go to the limit d —> 4 only at the end: 


^A{r) = f 


Of o 


1 + Cp— -6 log- 1-67 + 5 

47r V pT 


(14) 


where 7 is the Euler constant. Using the technique of [^|, it is not difficult 
to calculate the two-loop correction to However, we decided not to 

derive the next-to-next-to-leading term in the sum rule for 51 because 
the sum rule for is known only with the next-to-leading accuracy |]l3, 14]. 


As usual, we assume that the excited states’ contribution is dual to 


that of the theoretical spectral density (13) above the effective continuum 
threshold Ec- Equating the phenomenological expression for the correlator 
to the theoretical one at an imaginary time t = —i/E we obtain the sum 
rule 


,-.iE 

<qq> 5 mf rn^<qq> 

~ 3f^E 48/2^3 ’ 
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where 


X 


Em(x) 


/ 


1 — e * 


t 


dt 
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X X _|_ X^ 

m 2 ^ 3^ 

logx + 7 + 0{e~^/x), 


X <C 1 
X > 1 


(16) 


The /x-dependence of the right-hand side of (^) coincides with that of 
with the accepted accuracy. The correction appears to be about 
25 — 30% at ^ = 1 GeV and therfore is important. 



Figure 3: 5iF^/(380MeV)^ as a function of Borel parameter E: dashed 
line—without the perturbative correction, solid lines—with the perturbative 


correction (e = 500MeV Q, /j = IGeV, a* = 0.34, <qq> = —(240MeV)^, 
continuum threshold £c = IGeV, 1.25GeV, 1.5GeV). 

We can multiply (ll5| ) by exp(e/F) and differentiate in E in order to 
exclude c. This procedure is equivalent to finding a plato linear in the Borel 
parameter |17]. The result is shown in Fig. |^. The value of F is a subject 
of controversy in the current literature. We choose the middle value of |H] 

= (380MeV)^ as a reference point; the sum rule produces 51 F^, therefore 
lower values of F^ yield larger values of gi. The OFF series ( [l^ behaves 
as 1 + 270MeV/F — (170MeV/F)^ — (240MeV/F)^, so it seems that the 
expansion is applicable at F > 500MeV. This is close to the lower bound 


of the applicability region of the ordinary sum rule |12]. The continuum 
contribution is small until very large values of F, because the continuum 
only appears at the a^/vr level ([T^). At large F the sum rule ([T5|) (without 
the perturbative correction) leads to 


51 ^^ = ^<Qq>- 


(17) 









At lower E, stability of this sum rule is not very good. Our final numerical 
estimate is 5iF^/(380MeV)^ = 0.1 0.2. The effect of the perturbative 
correction is moderate, and dependence on the continuum threshold is very 
week. 

The numerical results of the previous works differ from each other largely 
due to different choices of the value of Expressed in a common scale, they 
are remarkably similar: 5iF^/(380MeV)^ = 0.16±0.03 |^], 0.16±0.04 p^] , 
0.17 ± 0.01 |5| (in the last result only diviation with Borel parameter is 
indicated and 20% level of accuracy is assumed), 0.14 ± 0.03 |@. The 


perturbative correction was not included in earlier works, and hence they 
had the continuum spectral density equal to zero. Ovchinnikov [^] included 
the first excited state in the phenomenological side, assuming its energy 
known. However, he only included 6{u; — e') in and not 6'{io — e'). 

We see no justification for such an assumption. Taking C 2 S'{uj — e') into 
account would add one more unknown parameter C 2 to the problem. 


3 Baryons 


Correlators of the SS and AS currents have the A-terms 

i<TjA{x)jhi'j'i0)>A = A^li,df,-^6{x)UA{xo). (18) 


If we define < 0 |j)^|A> = FacI, < 0 U 2 |S> = Tse^, then the physical 
states’ contributions to the correlators (Fig. ||) are 


nA(w) 

IIa{u}) 


I c I 

(a; — w — es 

2ff3FATs CA cs 

(cu - eA)(w - es) w-eA - et, 


(19) 


It is impossible to separate the 53 term from the mixed A-excited and S~ 
excited contributions unambiguously. We can do it approximately if A = 
es — ea ^ Gc — £a,e because in such a case partial fraction decomposition of 
the first term would give large contributions to ca,s ~ 1 /A with the opposite 
signs while the natural scale of ca,s in (H) is CA,s ~ l/(ec - eA,s)- This is 
not a defect of the sum rule but the uncertainty inherent to 53 which can 
be defined only when A is small compared to the chiral symmetry breaking 
scale. We choose to require ca = cs; the choices ca = 0 or cs = 0 would be 
equally good. 

The baryonic correlators in the OPE framework are described by the 
diagrams in Fig. The diagrams Fig. |^-c with the non-cut quark line 
interacting with the axial field vanish due to (P). We use the factorization 
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Figure 4: Correlator of baryon currents. 
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approximation for the four-quark condensate in Fig. In this approxi¬ 
mation two diagonal correlators <jiji> and <j 2 j 2 > coincide in both SS 
and AS cases, as well as <jh_\jT, 2 > and <jA 2 isi>- This is similar to the 


usual heavy baryon sum rules [lA, W], and confirms the observation that 
F\i = F \2 and = Fs 2 within the factorization approximation to the 
sum rules. Only even-dimensional condensates contribute to the diagonal 
SS and the nondiagonal AS correlators: 


nA(t) = 


2N\fl 


1 + 


5 CsX 'K‘^<qq>H‘^ 

3 ^ 6Nf^ 


( 20 ) 


where Cb/Cf = 1/(A^ — 1) (this term comes from the diagram Fig. ^). 
Only odd-dimensional condensates contribute to the nondiagonal SS and 
the diagonal AS correlators: 


nA(t) = 


2A^! <qq> 


2N 


1 - 


(l + IF 


2j.2 I 5 2^2\ 

mot +^m^t J 


( 21 ) 


These correlators correspond to the spectral densities 

N'.fl 


Pa{uj) = 


Ntt'^ 


^ ,5 Cb\ ml 

U ^ Cf) 6 


2N\<qq>u 


(plus (5(a;) and its derivatives). 

We use the standard continuum model = p^^^™'d{u}—ec)- Equat¬ 

ing the OPE (|^, ^) and the spectral representation aX t = —i/E, we obtain 
the sum rules 


( ‘^92 F^ 

V ^ 


+ c e 


-£e/E 


4f2 , 

_ J TT ^3 

^2 




X<qq> 




TT^ 


/i i^c/E) -|- 




I 1 _ "^0 _ \ 

2E2 I 16E2 36^2 ) 


(23) 


/ ‘^'PsFaEf, 

V A 

r2 


tanh — -h F 

2E 2 J ^ ' 


7m? 


= ^E^ h (sclE) - ^/o (ec/E) + 


7r2<gg>2 

18/2^4 


4:<qq> ^2 

2 ^ 
TT^ 




( 24 ) 


where fn{x) = 1 - e ^ Em=o x^lml. 
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Now we can start the numerical analysis. We adopt the energies and 
residues of A and S obtained from the HQET sum rules j 




eA = O.SGeV, F\ = 0.03GeV^, 
es = l.OGeV, Fs = 0.04GeV^ 

We assume that the effective continuum threshold for the A-terms of the S- 
S correlators is the same as without the external field: Ec = (1.3 ± 0.1)GeV; 
for the A“S correlators we use Ec = (1.1 ± 0.1)GeV. 

Fig. 1^ shows the results for g 2 F^/{l .6 ■ 10“^GeV®) and g^Ff^F^./{I-2 ■ 
10“^GeV®) obtained from the diagonal S-S and the nondiagonal A-S corre¬ 
lators (the first formulae in ( [^) and (^)). The spectral density (p^) behaves 
like and hence the continuum contribution to the sum rules ( P^ quickly 
grows: it is equal to the result at F ~ GOOMeV and is three times larger 
than the result at F ~ 900MeV. If we assume that the uncertainty in the 
continuum contribution is, e. g., 30%, then we can’t trust the sum rules at 
F > DOOMeV. The lower bound of the applicability region is determined 
by the convergence of the OPE series for the correlators. It behaves like 
1 — (190MeV/F)^ + (245MeV/F)^ for the S-S correlator; in the A-S case 
140MeV enters instead of 190MeV. It seems that OPE should be applicable 
at F > 400MeV. The sum rules in this window are unstable and therefore 
yield no definite predictions. 

Fig. |p3 shows the results for g 2 F^/(l.Q ■ 10“^GeV®) and g^Fi^F^,/{1.2 ■ 
10“^GeV®) obtained from the nondiagonal S-S and the diagonal A-S corre¬ 
lators (the second formulae in 


and (1^). The spectral density (^ be¬ 
haves like a;, and the continuum contribution grows not so quickly: it is equal 
to the result at F ~ IGeV and is three times larger than the result at F ~ 
l.GGeV. The OPE series behaves like l-|-(290MeV/F)^ [l — (150MeV/F)^], 
and the applicability region starts at a larger F. Stability of these sum rules 
is quite good. The sum rules give 52,3 = 0.4 0.7. 

In conclusion, the meson sum rule in an external axial field, with the 
perturbative correction and three successive power corrections, gives the 
value (7iF^/(380MeV)^ = O.l-i-0.2, which is much lower than the constituent 
quark model expectation gi = 0.75 unless a substantially lower value of F^ 
is used. The sum rules for g 2 ^'i following from nondiagonal S-S and diagonal 
Ar-Ti baryon correlators in an external axial field suggest 52,3 = 0.4 0.7, 

while diagonal S-S and nondiagonal A-S baryon sum rules have too large 
uncertainties. 
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